Generalized contact process with two symmetric absorbing states in two dimensions 
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We explore the two-dimensional generalized contact process with two absorbing states by means 
of large-scale Monte-Carlo simulations. In part of the phase diagram, an infinitesimal creation rate 
of active sites between inactive domains is sufficient to take the system from the inactive phase to 
the active phase. The system therefore displays two different nonequilibrium phase transitions. The 
critical behavior of the generic transition is compatible with the generalized voter (GV) universality 
class, implying that the symmetry-breaking and absorbing transitions coincide. In contrast, the 
transition at zero domain-boundary activation rate is not critical. 
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I. INTRODUCTION 

Phase transitions between different nonequilibrium 
steady states are a topic of great current interest in sta- 
tistical physics. These transitions display large-scale fluc- 
tuations and collective behavior over large distances and 
long times just as equilibrium phase transition. They oc- 
cur, for example, in surface growth, granular flow, chem- 
ical reactions, population dynamics, and even in traffic 
jams 0-01- 

The so-called absorbing state transitions are a partic- 
ularly well-studied type of nonequilibrium phase transi- 
tions. They separate fluctuating (active) steady states 
from absorbing (inactive) states where fluctuations stop 
completely. Generically, absorbing state transitions are 
in the directed percolation (DP) Q universality class; 
Janssen and Grassberger [t| [l(| conjectured that all ab- 
sorbing state transitions with a scalar order parameter 
and short-range interactions belong to this class as long 
as there are no extra symmetries or conservation laws. 
This conjecture has been confirmed in countless theo- 
retical and computer simulation studies. Experimental 
verifications were found in ferrofluidic spikes |11| and in 
the transition between two turbulent states in a liquid 
crystal fl2| . 

In recent years, significant attention has focused on ab- 
sorbing state transitions in universality classes different 
from DP that can occur if the system features additional 
symmetries or conservation laws. In 1997, Hinrichsen 
[la ] suggested several nonequilibrium stochastic lattice 
models with n > 2 absorbing states. In the case of two 
symmetric absorbing states (n = 2), he found the critical 
exponents to be different from the DP values. The corre- 
sponding universality class has been given several differ- 
ent names in the literature such as the ^-symmetric di- 
rected percolation class (DP2) or the directed Ising (DI) 
class. If the symmetry between the two absorbing states 
is broken, the critical behavior reverts back to DP. 

Recently, we revisited [3] one of the stochastic lat- 
tice models introduced in Ref. (la) , viz., the generalized 
contact process with two absorbing states in one space 
dimension. By employing large-scale Monte-Carlo sim- 
ulations, we found a rich phase diagram featuring two 



different nonequilibrium phase transitions separated by a 
special point that shares some characteristics with a mul- 
ticritical point. The generic transition occurs at nonzero 
values of the infection, healing and domain-boundary ac- 
tivation rates. It belongs to the above-mentioned DP2 
or DI universality class which in one dimension coincides 
[|| with the parity-conserving (PC) class [l5| (occurring, 
e.g., in the branching-annihilating random walk with an 
even number of offspring (BARWE) [16]). In addition, 
we found an unusual line of phase transitions at zero 
domain-boundary activation rate which turned out to be 
non-critical. 

Here, we consider the generalized contact process with 
two symmetric absorbing states in two space dimensions. 
The purpose of this paper is twofold. First, we wish 
to investigate whether the two-dimensional generalized 
contact process also displays the above-mentioned rich 
phase diagram having two nonequilibrium phase transi- 
tions. Second, we wish to study the critical behavior of 
these transitions and their universality. According to a 
conjecture by Dornic et al. [13], transitions with Zi sym- 
metry and no bulk fluctuations (i.e., transitions with two 
symmetric absorbing states) should be in the generalized 
voter (GV) universality class for which the upper criti- 
cal dimension is exactly two. Alternatively, the transition 
could split into a symmetry-breaking Ising transition and 
a DP transition [lj| Oil- To address these questions, we 
perform large-scale Monte-Carlo simulations. 

Our paper is organized as follows. We introduce the 
generalized contact process with several absorbing states 
in Sec. [II] Sec. Mil is devoted to the results and interpre- 
tation of our Monte-Carlo simulations. We conclude in 
Sec.llV] 



II. 



GENERALIZED CONTACT PROCESS WITH 
SEVERAL ABSORBING STATES 



We first define the simple contact process [20(, one 
of the prototypical models in the DP universality class. 
Each site r of a c?-dimensional hypercubic lattice can be in 
one of two states, either A, the active (infected) state or 
I, the inactive (healthy) state. During the time evolution 
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of the contact process, active sites infect their nearest 
neighbors, or they heal (become inactive) spontaneously. 
More rigorously, the contact process is a continuous-time 
Markov process during which active sites become inac- 
tive at a rate p, while inactive sites turn active at a rate 
Am/ (2c?) where m is the number of active nearest neigh- 
bor sites. The healing rate p and the infection rate A are 
external parameters. 

The long-time state of the contact process is deter- 
mined by the ratio of these two rates. If p ^> A, healing 
occurs much more often than infection. Thus, all in- 
fected sites will eventually become inactive, and the ab- 
sorbing state without any active sites is the only steady 
state. Consequently the system is in the inactive phase 
for /j > A. In the opposite limit, A 3> /i, the infection 
survives for infinite times, i.e., there is a steady state 
with a nonzero density of active sites. This is the active 
phase. These two phases are separated by a nonequilib- 
rium phase transition in the DP universality class occur- 
ring at some critical value of the ratio A///. 

Following Hinrichsen (l3| . we now generalize the con- 
tact process to n absorbing states. Each lattice site can 
now be in one of n + 1 states, the active state A or one 
of the n different inactive states [k = 1 . . . n). k is 
sometimes referred to as the "color" index. The Markov 
dynamics of the generalized contact process is defined via 
the following transition rates for pairs of nearest-neighbor 
sites, 

w(AA -» Alk) = w(AA ->■ I fe A) = p,/n , (1) 
w(AI k -> Ifelfe) = w(I k A -> 1*1*) = p k » (2) 
w(AI k -> AA) = w(I k A ->• AA) = A , (3) 
wihh -> I fc A) = 10(1*1, -► Ah) =a, (4) 

with k, I = 1 . . . n and k ^ I. All other transition rates 
vanish. We are mostly interested in the fully symmet- 
ric case, Hk = A* f° r a ^ For n = 1 and p, = p, the 
so defined generalized contact process coincides with the 
simple contact process discussed above. One of the rates 
/2, /i, A, and a can be set to unity without loss of gener- 
ality, thereby fixing the unit of time. We choose A = 1 
in the following. Moreover, to keep the parameter space 
manageable, we focus on the case p, = p [2l| . 

The rate ([4]) is responsible for the new physics in the 
generalized contact process. It prevents inactive domains 
of different color (different k) to stick together indefi- 
nitely. By creating active sites at the domain wall, the 
two domains can separate. Thus, the rate §4§ allows the 
domain walls to move through space. We emphasize that 
without the process i.e., for a = 0, the color of the 
inactive sites becomes unimportant, and all Ifc can be 
identified. Consequently, for a — 0, the dynamics of the 
generalized contact process reduces to that of the simple 
contact process for all values of n. In the main part of 
this paper, we shall focus on the case of n = 2 inactive 
states. 

Before we turn to our Monte-Carlo simulations of 
the two-dimensional generalized contact process, let us 



briefly summarize the simulation results in one dimen- 
sion 14] for comparison. For a = 0, i.e., in the absence 
of the boundary activation process ((4]), the system un- 
dergoes an absorbing state transition at a healing rate 
/i = p c c p rs 0.303, which agrees with the critical heal- 
ing rate of the simple contact process. In agreement 
with the general arguments above, this transition is in 
the DP universality class. For healing rates between p cp 
and p* » 0.552, the system is inactive if a — but an 
infinitesimal nonzero a takes it to the active phase. Fi- 
nally, for p > p* , the transition occurs at a finite nonzero 
value of a. The one-dimensional generalized contact pro- 
cess with two inactive states thus has two lines of phase 
transitions, (i) the generic transition occurring at p > p* 
and a = u c {p) > and (ii) the transition occurring for 
p cp < p < p* as a approaches zero. 



III. MONTE-CARLO SIMULATIONS 
A. Method and phase diagram 

In order to address the two main problems raised in 
the introduction, viz, the phase diagram of the two- 
dimensional generalized contact process with two inactive 
states and the critical behavior of its phase transitions, 
we performed two types of large-scale Monte Carlo sim- 
ulations, (i) decay runs and (ii) spreading runs. Decay 
runs start from a completely active lattice; we measure 
the time evolution of the density p(t) of active sites as 
well as the densities pi(t) and p2(t) of sites in inactive 
states Ii and I2 , respectively. Spreading simulations start 
from a single active (seed) site embedded in a system of 
sites in state Ii . Here we monitor the survival probability 
P s (t), the number of sites in the active cloud N s (t) and 
the mean-square radius of this cloud, R 2 (t) . 

In both types of runs, the simulation is a sequence 
of individual events. In each event, a pair of nearest- 
neighbor sites is randomly selected from the active region. 
For the spreading simulations, the active region initially 
consists of the seed site and its neighbors; it is updated 
in the course of the simulation according to the actual 
size of the active cluster. For the decay runs, the active 
region comprises the entire sample. The selected pair 
than undergoes one of the possible transitions according 
to eqs. (HJ) to (fj} with probability tw. Here the time step 
r is a constant which we fix at 1/2. The time increment 
associated with the event is r /N pa i r where N pa i r is the 
number of nearest-neighbor pairs in the active region. 

Using this procedure, we investigated the parameter 
region 0.5 < p < 1.2 and < a < 1. We simulated 
samples with sizes up to 20000 x 20000 sites for times 
up to t rnax = 3 x 10 6 . The a — p phase diagram that 
emerged from these calculations is shown in in Fig. [TJ 

In many respects, it is similar to the phase dia gra m of 
the one-dimensional generalized contact process [14j . In 
the absence of the domain-boundary activation process 
(i.e., for er = 0), the transition from the active phase to 
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FIG. 1. (Color online) Phase diagram of the two-dimensional 
generalized contact process with two inactive states as func- 
tion of the healing rate fi and the domain-boundary activation 
rate a. For [i < ^ — 0.6066, the system is in the active phase 
for any a. For < [i < (i* = 1.0000, the system is inactive 
at a = (thick solid red line), but an infinitesimal a takes 
it to the active phase. For /i > fi* , the system is inactive for 
any a . 



the inactive phase occurs at a healing rate of \i = \xf = 
0.6066(2) which agrees well with the critical point of the 
simple contact process (see, e.g., Refs. [IH HH). For heal- 
ing rates in the interval [if < fi < [l* = 1.0000(2), the 
generalized contact process is inactive at a — 0, but an 
infinitesimal nonzero a takes it to the active phase. Thus, 
we find a line a phase transitions at fif < \i < fi* and 
(7 = 0. In addition to this line of a — absorbing state 
transitions, we also find a line of generic (nonzero a and 
fj,) transitions. In contrast to one space dimension, this 
line is exactly "vertical" within our accuracy, i.e., the 
critical healing rate fi c = 1.0000(2) does not depend on 
a for all a > 0. We note in passing that our critical heal- 
ing rate is in agreement with the estimate fi c « 0.99(1) 
obtained in Ref. [HI for a = 1. 

In the following subsections we shall discuss in detail 
the properties of both phase transition lines as well as 
special point {[i* , 0) that separates them. 



FIG. 2. (Color online) Spreading simulations at a = 0.1 for 
several /i close to the phase boundary. Main panel: Survival 
probability P 3 as a function of time t. Inset: Number N a of 
active sites as a function of time t. The data close to criticality 
are averages over 10 6 runs on a 4000 x 4000 system, smaller 
numbers of runs were used away from criticality. 
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B. Generic transition 

In order to identify the generic transition and to study 
its critical behavior, we performed sets of spreading sim- 
ulations at constant domain-boundary activation rate 
a = 0.01, 0.05, 0.1, 0.5 and 1. For each a, we have varied 
the healing rate \i from 0.8 to 1.1. Figure [5] shows the re- 
sulting time evolution of the survival probability P s and 
the number of sites in the active cloud N s (t) for a = 0.1 
and several fi. The data indicate a critical healing rate of 
/j, c = 1.0000(2) for this a value. Analogous simulations 
for a — 0.01,0.05,0.5 and 1 yielded, somewhat surpris- 
ingly, exactly the same critical healing rate. We thus 



FIG. 3. (Color online) Survival probability P a and number of 
active sites N 3 as functions of t for several points located on 
the generic phase boundary /i = 1.0000 (2 x 10 6 to 10 7 runs 
used). Inset: prefactor B„ vs. a. The straight line is a fit to 
a power-law B a ~ cr - ''. 



conclude that in the two-dimensional generalized contact 
process, the critical healing rate [l c is independent of a 
for all er > 0. 

Figure [3] shows the survival probability P s and number 
N s of active sites as functions of time for all the respective 
critical points. In log-log representation, the long-time 
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FIG. 4. (Color online) Survival probability P s (t) for several 
points located on the generic phase boundary plotted as t P s 
vs. ln(t). Straight lines correspond to mean-field behavior 
with logarithmic corrections. Inset: Same data plotted as 
\n(tP s ) vs. ln(t). Straight lines represent pure power laws. 



FIG. 5. (Color online) Density of active sites plotted as p - 
vs. ln(t) for several points located on the generic phase bound- 
ary. The data are averages over 100 runs with system size 500 
x 500. The curve for a = 0.01 is shown in the inset because 
its density values are much smaller than those of the other 
curves. 



parts of the N s and P s curves for different a are perfectly 
parallel within their statistical errors, i.e., they differ only 
by constant factors, confirming that the critical behavior 
of the generic transition is universal. Fits of the long- 
time behavior to the pure power laws P s = B a t~ s and 
N s = C a t e give estimates of S = 0.900(15) and 6 = 
—0.100(25). These values are very close to the mean- 
field values 8m f = 1 and @mf — 0. According to the 
conjecture by Dornic et al. [17], the generic transition 
should be in the GV universality class. Because the upper 
critical dimension of this universality class is exactly two, 
this conjecture corresponds to mean-field behavior with 
logarithmic corrections. 

To test this prediction we compare in Fig. 2] plots 
of ln(iP s ) vs. ln(t) (straight lines corresponds to power 
laws) and t P s vs. ln(i) (straight lines correspond to log- 
arithmic behavior). Although both functional forms de- 
scribe the long-time data reasonably well, the curves in 
the \n(tP s ) vs. ln(i) plot show a systematic downward 
curvature. Moreover, the semi-logarithmic plot, tP s vs. 
ln(t), leads to straight lines over a longer time interval 
which we take as evidence for GV critical behavior. We 
performed an analogous analysis for number of active 
sites N s . Again, both a simple power law and mean- field 
behavior with logarithmic corrections describe the data 
reasonably well, with the quality of fits being somewhat 
higher for the latter case. We also measured (not shown) 
the mean-square radius R 2 {t) of the active cloud as a 
function of time. A pure power-law fit of its long time 
behavior, R 2 (t) ~ t 2 '\ gives 2/z = 0.97(4) (z = 2.06(8)). 
The data can be described equally well by mean-field be- 
havior R 2 (t) ~ t with logarithmic corrections. 

In addition to the spreading runs, we also performed 
density decay runs at the generic phase boundary. The 
resulting density of active sites p as a function of time 
can be fitted with a pure power law p(t) ~ t~ a giv- 



ing a very small value of a — 0.080(4). A better fit is 
achieved with the simple logarithmic time dependence 
p(t) r~j l/ln(t/to) (with to a microscopic time scale) ex- 
pected for the GV universality class. This type of behav- 
ior is demonstrated in Fig. [5j 

In summary, although all our results for the generic 
transition can be fitted both by pure power laws and by 
mean-field behavior with logarithmic corrections, the lat- 
ter functional forms yield fits of somewhat higher quality. 
We also note that the critical exponents resulting from 
the pure power-law fits approximately fulfill the hyper- 
scaling relation <d — d/z = —a — S. However, the agree- 
ment is not very good (in particular, it is significantly 
worse than in one dimension [Hj]), indicating that the 
measured pure power-laws are not the true asymptotic 
behavior. Our results thus support the conjecture that 
the generic transition of the two-dimensional generalized 
contact process with two inactive states is in the GV uni- 
versality class. 

C. Transition at a — 

After addressing the generic transition, we now discuss 
in more detail the line of phase transitions occurring at 
cr = and pf < p < p* . To study these transitions, we 
carried out several sets of simulations for fixed healing 
rate p and several a values approaching a = 0. 

We start by discussing the density decay runs. Figure 
[5] shows the stationary density p s t of active sites (reached 
at long times) as function of a for several values of the 
healing rate p. The figure shows that the stationary den- 
sity depends linearly on a for all healing rates in the in- 
terval p^P < p < p* , i.e., p st — B^a^ with u = 1 and 
i? M being a ^-dependent constant. We also analyzed how 
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FIG. 6. (Color online) Density decay simulations. Main 
panel: stationary density p 3t as a function of the boundary 
rate a for various healing rates p. For p^, p </*</**, the 
solid lines are fits of the low-a behavior to p s t = i3 M <r. At 
the simple contact process critical point, p — p cp — 0.6066, 
we fit to the power-law p 3t ~ u" cp which gives an exponent 
of Lu C p — 0.274(5). The data are averages over 300 to 600 
runs with system sizes 100 x 100. Inset a: prefactor B M of 
the linear a dependence as a function of p — p cp . A fit to a 
power law gives ~ (p — Pc P )~ K with k = 1.56(5). Inset 
b: prefactor as a function of p* — p. A fit to a power law 
gives B^ ~ (fi* — p) K with k* ~ 0.23. 

the prefactor of this mean-field-like behavior depends 
on the distances from the simple contact process critical 
point and from the special point at p — ji* and a = 0. 
As inset (a) of Fig. [6] shows, diverges as (/z — Pc P )~ K 
with k = 1.56(5). According to inset (b), it vanishes as 
(/i* — p) K with k* w 0.23 when approaching p* . 

At the critical healing rate p\ cp of the simple con- 
tact process, the stationary density displays a weaker 
tr-dependence. A fit to a power-law p st ~ tr" cp gives 
an exponent value of uj cp = 0.274(5). 

Let us now compare these results with the behavior 
of spreading simulations in the same parameter region. 
Figure [Jj shows the survival probability P s (t) and the 
number of active sites N s (t) for a fixed healing rate of 
fi = 0.8 and several values of the boundary rate a. After 
an initial decay, the number of active sites grows with 
time for all a values, establishing that the system is in 
the active phase for all a > 0. In agreement with this, 
the survival probability approaches a nonzero constant 
in the long-time limit. Remarkably, this stationary sur- 
vival probability does not approach zero with vanishing 
a. Instead, it approaches a cr-independent constant. We 
performed similar sets of simulations at other values of /j, 
in the range ix cp < /i < fi* , with analogous results. 
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FIG. 7. (Color online) Spreading simulations: Survival prob- 
ability P„ and number of active sites N 3 as functions of time 
t for a fixed healing rate of /i = 0.8 and several a. The data 
are averages over 2000 to 10000 runs on a 4000 x 4000 system. 



We thus conclude that the behavior at the a — transi- 
tion of the two-dimensional generalized contact process is 
very similar to the one-dimensional case. It can be under- 
stood in terms of the domain- wall motion as follows [3] . 
The relevant long-time degrees of freedom at p, > [j%P and 
u«l are the domain walls between Ii and I2 domains. 
These walls can hop, branch and annihilate. The crucial 
observation is that the rates which control the domain 
wall dynamics are all proportional to a for a <C 1, im- 
plying that their ratios are cr-independent. Consequently, 
the stationary state of the domain walls does not depend 
on a for a <C 1. This explains why the survival proba- 
bility P s saturates at a nonzero, u-independent value in 
Fig. [7] It also explains the c-dependence of the station- 
ary density p s t because active sites are created mostly at 
the domain walls at rate a. Therefore, their stationary 
density is proportional to both a and the stationary do- 
main wall density pdw, i-e., p s t ~ o~Pdw, in agreement with 
Fig. [SJ Based on this argument, the exponent k* in inset 
(b) of Fig. [5] should be identical to the exponent /3 of the 
generic transition line flij ]. which vanishes in mean- field 
theory. Our value, k* w 0.23 is thus somewhat too high 
which we attribute to it not representing the asymptotic 
behavior, in agreement with the significant curvature of 
the data in inset (b) of Fig. [6] 

Just as in one dimension, the phase transition line at 
(7 = and pic P < < M* is thus not a true critical line. It 
only appears critical because the stationary density p s t 
(trivially) vanishes with a. Correspondingly, the time 
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evolution right on the transition line a = does not 
display critical power laws. This also implies that the 
point (/J,, a) = (p*,0) is not a multicritical point, but a 
simple critical point in the same universality class as the 
generic transition. 

D. Scaling of p st at the contact process critical 
point (p c c p ,0) 

The behavior of the stationary density of active sites 
p s t close to the simple contact process critical point at 
p = pf and a — can be understood in terms of a 
phenomenological scaling theory. We assume the homo- 
geneity relation 

p st {Ap, er) = b^Kp^Apb^K^b^) (5) 

where Ap = p, — p cp , and 6 is an arbitrary scale fac- 
tor. f3 cp = 0.584 and v^ p = 0.734 are the usual order 
parameter and correlation length exponents of the two- 
dimensional contact process [22|, [23[, and y cp denotes 
the scale dimension of er at this critical point. Setting 
b = erVacp gives rise to the scaling form 

p st (Ap 7 a) = <j^/ { < v ^X (Ap<j- l/{u ^ v ^ (6) 

where X is a scaling function. At criticality, Ap = 0, this 
leads to p st (0,a) - a ^ P /(^ p Vc P ) ( using X (0) = const). 
Thus, w cp = (3 cp / '(v^pijcp). For a — > at nonzero Ap, 
we need the large-argument limit of the scaling function 
X. On the active side of the critical point, Ap < 0, the 
scaling function behaves as X{x) ~ \x\^ cp to reproduce 
the correct critical behavior of the density, p st ~ \p — 

H-c \ 

On the inactive side of the critical point, i.e., for 
Ap > and a — > 0, we assume the scaling function to 
behave as X(x) ~ x~ K '. We thus obtain p st ~ (Ap,)~ K o- u 
(just as observed in Fig. [6]) with ui — (f3 cp + k) / (v^r p y cp ). 
As a result of our scaling theory, the exponents lo,uj cp 
and k are not independent, they need to fulfill the rela- 
tion w C p(/3 cp + k) = PcpLJ. Our numerical values, u = 1, 
u C p — 0.274 and n — 1.56 fulfill this relation in very good 
approximation, indicating that they represent asymp- 
totic exponents and validating the homogeneity relation 
([5]) . The resulting value for the scale dimension y cp of a at 
the simple contact process critical point is y cp = 2.9(1). 

IV. CONCLUSIONS 

To summarize, we investigated the two-dimensional 
generalized contact process with two inactive states by 
means of large-scale Monte-Carlo simulations. Its global 
phase diagram is very similar to that of the correspond- 
ing one-dimensional model. In particular, the generic 
(a > 0) phase boundary between the active and inac- 
tive phases does not continuously connect to the criti- 
cal point of the a = problem, i.e., the critical point 



(Pc P , 0) of the simple contact process. Instead, it termi- 
nates at a separate end point {p*,0) on the p, axis. As 
a result, the two-dimensional generalized contact process 
has two nonequilibrium phase transitions. In addition to 
the generic transition occurring for a > 0, there is a line 
of transitions at a = and pf p < p < p* . We note that 
there is one interesting difference between the phase dia- 
grams in one and two dimensions. In one dimension, the 
critical healing rate p, c increases with increasing bound- 
ary rate a. In contrast, the results of this paper show 
that the critical healing rate in two dimensions is com- 
pletely independent of a. Moreover, its value seems to 
be equal to unity (i.e., equal to that of the infection rate 
A). The reason for this peculiar behavior is presently an 
open question. 

To determine the critical behavior of the generic tran- 
sition, we performed simulations at and close to several 
points on the generic (a > 0) phase boundary. We found 
the same critical behavior for all of these points, i.e, it 
is universal. Our data can be fitted reasonably well with 
pure power laws, giving the exponents O = —0.100(25), 
S = 0.900(15), a = 0.080(4), and z = 2.06(8). How- 
ever, fits of equal and sometimes even better quality over 
longer ranges of time can be obtained by fitting to mean- 
field critical behavior, O = 0, 6 = 1, a = 0, and z = 2 
with logarithmic corrections. Our results thus support 
the conjecture [l7[ that the critical behavior of the two- 
dimensional generalized contact process is right at its 
upper critical dimensions. (This implies that the DP2 
class in two dimensions coincides with the GV class.). 
We also note that our simulations showed no indications 
of the transition being split into a symmetry-breaking 
transition and a separate DP transition as found in some 
absorbing-state Potts models [l8| . 

As in one space dimension, the line of transitions at 
(j = and Pc P < p < p* is not a critical line. The survival 
probability P s remains finite when approaching this line. 
The density p of active sites vanishes, but simply because 
the domain-boundary activation rate a vanishes. The 
behavior in the vicinity of the transition line is controlled 
by the dynamics of the I1-I2 domain walls which is not 
critical for pf < p < p* . 

Crossovers between various universality classes of ab- 
sorbing state transitions in one dimension have been in- 
vestigated by several authors [24T - [27l | . Some of the sce- 
narios lead to conventional crossover scaling (of the type 
ct c ~ (p- /^ p ) 1/0 ). Park and Park [H[ found a discon- 
tinuous jump in the phase boundary along the so-called 
excitatory route from infinitely many absorbing states to 
a single absorbing state. There also is some similarity 
between our mechanism and the so-called channel route 
26] from the PC universality class to the DP class which 
involves an infinite number of absorbing states charac- 
terized by an auxiliary density (which is density of I1-I2 
domain walls in one-dimensional the generalized contact 
process [Hj])- To the best of our knowledge, a similarly 
systematic investigation of crossovers between absorbing 
state universality classes in two space dimensions has not 



yet been performed. 

As our results suggest that the two-dimensional gener- 
alized contact process is right at the upper critical dimen- 
sions, the critical behavior of its (generic) phase transi- 
tion in dimensions d > 2 should be governed by mean- 
field theory. 
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